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Abstract 

We present a systematic hydrodynamic study of the evolution of hadron spectra and their az- 
imuthal anisotropy from the lowest collision energy studied at the Relativistic Heavy Ion Collider 
(RHIC), V* = 7.7 A GeV, to the highest energy reachable at the Large Hadron Collider (LHC), 
V* = 5500 A GeV [4J. The energy dependence of the flow observables are quantitatively studied 
for both the Monte-Carlo Glauber and Monte-Carlo Kharzeev-Levin-Nardi (MC-KLN) models. 
For MC-Glauber model initial conditions with rj/s = 0.08, the differential charged hadron elliptic 
flow v^ipj, ~\fs) is found to exhibit a very broad maximum in the region 39 < < 2760 A GeV. 
For MC-KLN initial conditions with rj/s = 0.2, a similar "saturation" is not observed up to LHC 
energies. We emphasize that this "saturation" of elliptic flow arises from the interplay between 
radial flow and elliptic flow which shifts with tJs depending on the fluid's viscosity. By gener- 
alizing the definition of spatial eccentricity to isothermal hyper-surface, we also calculate e x on 
the kinetic freeze-out surface at different collision energies. 



1. Introduction 

The recent Beam Energy Scan (BES) program HIEIO at the Relativistic Heavy Ion Collider 
(RHIC) pursues one of the major goals of heavy-ion collision experiments: to explore the QCD 
phase diagram and search for the phase boundary between the normal nuclear matter and quark- 
gluon plasma (QGP). The BES program at RHIC together with Pb+Pb collisions at the Large 
Hadron Collider (LHC) provide us with a unique opportunity to study systematically the collision 
energy dependence of a large number of relativistic heavy-ion collision observables. 

Here we study the collision energy dependence of charged hadron transverse momentum 
spectra and elliptic flow coefficients [4], using (2+l)-d viscous hydrodynamics coupled with a 
modern lattice QCD based equation of state |5]|6|. Our work focuses on qualitative tendencies 
rather than quantitative comparison with experimental data. It makes simplifying assumptions 
that are justified at high energies but gradually break down at lower longitudinal boost invari- 
ance, an equation of state for matter with zero net baryon density, and a purely hydrodynamic 
approach with constant specific shear viscosity, including for the dilute hadronic rescattering 
stage which at lower energies occupies an increasing fraction of the fireball's dynamical history 
and should be described microscopically [7 1. These limitations can be cured in future work; they 
are expected to modify our conclusions quantitatively but not qualitatively. 

2. Results and discussion 



Evolution of charged hadron multiplicity and total elliptic flow: In Figs.[TJa,b) we show the 
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Figure 1: (a): Cen- 
trality dependence of 
final charged hadron 
multiplicity per par- 
ticipant nucleon pair 
as a function of N vM 
for MC-Glauber initial 
conditions, with colli- 
sion energies varying 
from VJ = 7.7AGeV to 
VJ = 2760AGeV. (c): 
Centrality dependence 
of -jg* from the 
lower energy runs in 
(a) scaled up to the 
LHC results, for shape 
comparison. (b, d): 
Same as (a, c) but 
for MC-KLN initial 



200 



300 



200 



300 



400 conditions. 



centrality dependence of the charged hadron multiplicity for both MC-Glauber and MC-KLN 
models with collision energies from V* = 7.7 to 2760 A GeV. The curves account for viscous 
entropy production during the hydrodynamic evolution. At LHC and top RHIC energies our 
results for both initialization models agree well with the experimental data 18,9]. Our lower col- 
lision energy predictions can be checked against data from the RHIC BES program. In Figs.[TJ;,d 
we scale the lower energy results by constant factors to align them with the LHC curve in central 
(0-10%) collisions, to see how the centrality dependence changes with -\fs. For the MC-Glauber 
model, the curves fall almost on top of each other. This is because we keep the mixing ratio be- 
tween the wounded nucleons and binary collisions fixed for the low energy runs at RHIC, and it 
also reflects the fact that viscous entropy production is small and has little effect on the centrality 
dependence. For the MC-KLN model, however, the slope of the centrality dependence gets flat- 
ter as the collision energy decreases. Only the top RHIC and LHC energy curves approximately 
fall on top of each other. We found that this tendency originates in the nature of the MC-KLN 
model itself. Our MC-KLN calculations thus predict a violation of the y^-scaling of the central- 
ity dependence of at lower collision energies that is not seen with the MC-Glauber initial 
conditions. This may help to discriminate experimentally between these models. 

In Fig. [2] we explore the scaling of elliptic flow with charged hadron multiplicity density 
("multiplicity scaling") over a wider range of V* than previously studied, for both of the ini- 
tialization models. For MC-Glauber initial conditions (Fig. [2^) the "multiplicity scaling" curve 
V2/S2 vs. {l/S){dN c h/drj) shows excellent universality over the entire collision energy range be- 
tween 7.7 and 2760 A GeV. But for MC-KLN (Fig. [2J5), lower collision energies result in larger 
V2/f2 values at the same charged hadron multiplicity density (as previously shown in |10|). We 
find that the main reason for the different collision energy dependences of the two models lies in 
their different behavior of the initial overlap area S : As the collisions become more peripheral, S 
decreases more rapidly in the MC-KLN model than in the MC-Glauber model [4|. This slightly 
faster drop of S in the MC-KLN model shifts the "universal" scaling curves in Fig. [2] to the 
right and shrinks the covered range in (l/S)dN c h/dr]. The different V^-dependences of V2/f2 as 
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Figure 2: Eccentricity- 
scaled pr-integrated V2 
plotted as a function 
of the charged hadron 
multiplicity density for 
different collision ener- 
gies, for MC-Glauber 
initial conditions with 
ri/s = 0.08 (a) and 
MC-KLN profiles with 
rj/s = 0.2 (b), respec- 
tively. 
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Figure 3: (a, b): 
Differential elliptic flow 
of all charged hadrons 
at 20-30% centrality 
in Au+Au and Pb+Pb 
collisions for different 
collision energies, (c, 
d): yi-dependence of 
the differential charged 
hadron elliptic flow 
V? (j>t, V*) at 5 fixed 



Pt values 
GeV/c. 
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function of dN^/drj in Figs. [2^ and[2j3 thus reflect primarily the fact that the shape of the initial 
profiles evolves differently with centrality in the two initialization models. 

"Saturation" of differential elliptic flow: In Figs.[3^,b, we show the differential charged 
hadron elliptic flow for Au-Au or Pb-Pb collisions at 20-30% centrality. The differential el- 
liptic flow is affected by both total elliptic and radial flow. With MC-Glauber initial conditions 
the differential elliptic flow for pj < 2 GeV/c remains almost unchanged for ^[s > 39 A GeV. As 
the collision energy increases, both radial and elliptic flow increase, due to the longer fireball life- 
time. Naively the increasing total elliptic flow should also lead to a larger differential Vj, but this 
tendency is counteracted by the growing radial flow which blueshifts the momentum anisotropy 
to larger pj. For the runs with MC-KLN initial conditions we use a larger r\j s value. The result- 
ing larger viscous effects suppress the total elliptic flow at lower collision energies more strongly 
than for the MC-Glauber runs, leading to a monotonous decrease of the slope of the differential 
v i(Pt) with decreasing collision energy. To further illustrate this point we plot in Fig. |3};,d the 
Vs-dependence of v^ipr) at 5 fixed pj points. In this representation one sees that for the MC- 
Glauber model with 77/.? = 0.08, vf 1 at any fixed pj <2GeV/c features, as a function of s/s, a 
very broad maximum around top RHIC energy (200 A GeV). For lower p T < 0.5 GeV/c, heavier 
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Figure 4: Vs-dependence of the final spatial 
eccentricity £f of the isothermal kinetic freeze- 
out surface at T iec = 120 MeV, for 10-30% cen- 
trality. The initial eccentricity is 0.26 for the 
MC-Glauber model and 0.32 for the MC-KLN 
model. The experimental points indicate prelim- 
inary data 111! from an azimuthal HBT analysis 
by the STAR Collaboration. 



particles, or smaller tj/s this maximum moves towards lower V^- With the larger rj/s = 0.2 used 
in the MC-KLN model, the strong reduction of vS h at low collision energies shifts the maximum 
of vf at any fixed p T towards higher yfs; Fig. 3d shows that for rj/s = 0.2 this observable does 
not peak below the top LHC energy, except for very small pj < 200 Me V/c. 

Freeze-out shape analysis: In Fig. [4] we show the final fireball eccentricity calculated along 
the kinetic freeze-out surface, - 120 MeV, as a function of collision energy. As the collision 
energy increases, the final spatial eccentricity Sf decreases monotonically for both MC-Glauber 
and MC-KLN models. This is because at higher collision energy the fireball has longer lifetime 
to decompress its original deformation and become more isotropic. In Fig. [4] we compare our re- 
sults with recent STAR data from an azimuthal HBT analysis 1 1 1 1. Our calculations qualitatively 
agree with the experimental data and reproduce the trend of the collision energy dependence of 
Sf. MC-Glauber runs with tj/s — 0.08 quantitatively reproduce the data at -\fs = 200 A GeV while 
underpredicting the final eccentricity by ~10% at lower energies. MC-KLN initial conditions 
with tj/s — 0.2 result in 15-20% larger final eccentricities, due to the ~20% larger initial eccen- 
tricities of the MC-KLN profiles. Within the explored range, we found weak sensitivity of these 
curves to 77/ s. Extending our calculations to LHC energy we predict that Sf will approach zero 
around yfs = 2.76-5.5 A TeV 
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